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Abstract. The radio emission of a galaxy consists of 
thermal bremsstrahlung, synchrotron emission from dis- 
crete supernova remnants, and diffuse synchrotron emis- 
sion from cosmic ray electrons spread over the galactic 
disk and halo. Each of these components has a different 
spectral index so that the total radio spectral index of a 
galaxy depends sensitively on their relative contribution 
and on the processes shaping the diffuse synchrotron emis- 
sion. In the present paper we calculate the contribution of 
supernova remnants to the total synchrotron emission and 
conclude that it is about 10%. This moderate contribution 
has a noticeable effect on the nonthermal spectral index, 
lowering it by ~ 0.1 for steep spectra. We then calculate 
the diffuse synchrotron emission in two simple models, a 
diffusion and a convection model. We find that the spa- 
tially integrated nonthermal spectral index is in general a 
poor diagnostic for the type of propagation or the impor- 
tance of energy losses. Spatially resolved radio data for 
the halo of galaxies are necessary in order to draw firm 
conclusions. The steepening of the spectrum away from 
the disk is a clear indication that synchrotron and inverse 
Compton losses are taking place during the propagation 
of cosmic ray electrons in the halo. In those galaxies for 
which spatially resolved data for the halos exist such a 
steepening has been found. We conclude therefore that 
energy losses are generally important and that cosmic ray 
electrons cannot freely escape from galaxies. 
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1. Introduction 

The total radio emission of a galaxy consists of different 
components: The thermal radio emission, Pthorm, repre- 
senting bremsstrahlung from H II regions and the nonther- 
mal radio emission which is synchrotron emission from 
cosmic ray electrons (CREs). As far as the nonthermal 
radio emission, P n th, is concerned, a further division can 



be made by distinguishing between the diffuse radio emis- 
sion, -Pdiff, from CREs that are spread over the galactic 
disks and halos and the emission from discrete supernova 
remnants (SNRs), Psnr- The three components possess 
very different spectral indices: Whereas the thermal radio 
emission has a flat spectral index (oc ^^ 01 ), the nonther- 
mal radio spectrum is steeper. The spectral index of SNRs 
is on average w v~°- 5 (Green 1998) and the diffuse radio 
emission can have, depending on the dominant process 
for the CR propagation, spectral indices^ , a, between 
about 0.5 and 1.1 which can, additionally, vary over the 
frequency range. Therefore, the total radio spectral index 
depends very sensitively on the relative contribution of the 
different components. 

The spectrum of the diffuse synchrotron emission is 
shaped by the physical processes that are characterizing 
the CR propagation, in particular the type of propagation 
(diffusion or convection) , of energy losses (synchrotron, in- 
verse Compton, adiabatic losses or bremsstrahlung), and 
the confinement of CREs. Here, we consider CREs con- 
fined to a galaxy if they loose their energy through com- 
bined synchrotron and inverse Compton losses, down be- 
low the energy level corresponding to the observing fre- 
quency, before they are able to leave the galactic halo. The 
question whether CREs are confined to galaxies or can es- 
cape from the halos without suffering considerable energy 
losses is important for the interpretation of the far-infrared 
(FIR)/radio correlation. Different models have been pro- 
posed, explaining this correlation by advocating either the 
sit uation that CREs are confined (Volk [1989 ; Lisenfeld et 
al. 1996a ) or t hat th ey can escape m ore or less freely (Chi 

Bicay |1993| ; Niklas & Beck 



& Wolfcndalc 1990 Hclou 



19971) 



The interpretation of the radio spectrum in terms of 
the above mentioned processes is difficult because they 
are hard to disentangle without spatially resolved obser- 
vations. Furthermore, it is necessary to separate the dif- 
fuse radio emission from the thermal radio emission and 
the radio emission from SNRs. The observational sepa- 
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ration of the three contributions based on their different 
spectral indices is difficult, especially for the radio emis- 
sion of SNRs which has a similar spectral index as the 
diffuse synchrotron emission. The separation of the ther- 
mal radio emission is possible, if high- frequency radio data 
(y « 10 — 20 GHz) are available, because its spectral index 
is very different. 

In this paper we are going to examine the influence of 
these various processes on the radio spectral index. For 
this, we will start by estimating the contribution of SNRs 
to the radio emission and to calculate its effect on the 
spectral index. Then, we discuss within a simple model 
the influence of energy losses, the type of propagation and 
escape of CREs on the spectral index. 

2. Radio emission of SNRs 

2.1. Theoretically expected value 

We estimate the fraction of the nonthermal radio emission 
of a galaxy that is expected to come from SNRs in the 
framework of a simple leaky-box model. The number of 
CREs in a certain energy range dE, N(E), can then be 
calculated as: 



dN(E) N(E) 
dt + t 



Q(E). 



(1) 



Q(E) is the production rate of CREs and r the life-time 
of CREs. The main sources of CREs are gravit ational su- 
pernovae (SN) (SN II and Ib /c) (e.g. Xu et al. |l994a| and 
discussion in Xu et al 1994b ). The energy spectrum of the 
CREs produced can be described by a power law: 



(){ E) = f'sN '/s\ ( -Jj- 



(2) 



where i/sN is the SN-rate, <7sn the number of CREs pro- 
duced per energy interval per SN and 7 the energy in- 
jection spectral index. For CREs in SNRs r is given by 
the life-time of the SNR, tsnr, and for the diffuse CREs 
t = T diS can be expressed as: 



1 

Tdiff 



1 



1 



Tlo 



(3) 



where r osc is the escape time-scale from the galaxy and 
Tioss is the energy loss time-scale. We take into account a 
possible energy dependence of Tdiff according to: 



E x ~ v 
Tdi&(E) = r diff (E a ) ( — 



(4) 



The exponent y can have different values depending on 
the characteristics of the propagation of CREs. It is y = 
1 if Tdiff = T"ioss an d synchrotron and inverse Compton 
losses are the dominant energy loss processes. If Tdiff = 
t csc and CREs propagate mainly by diffusion, the energy 
dependence of the diffusion coefficient, D(E) — D E~^, 



gets reflected in y — fi ss 0.5. If, on the other hand, CREs 
can escape freely (Tdiff = T esc ) by convection or if adiabatic 
and bremsstrahlung losses are dominant, we expect y = 0. 

In the case of a steady state (ON /dt = 0) the solution 
of eq. (Ill) is simply: 



N(E) = Q(E) t. 



(5) 



We make the simplification that an electron with en- 
ergy E emits all the synchrotron radiation at the fre- 
quency 



0.435 



E 
m e c 2 



v G = 4.65 



f_^WA N )MHz (6) 
VlGeVy \(iGj 



(Longair 1992, Vol. 1 p. 332) where vq, is the gyrofre- 
quency. Then the total synchrotron emission is given by: 



P{v)dv = N(E) 



dE 



dE 



dv, 



(7) 



with the energy losses by synchrotron radiation 



dE 

I = — (TtC 



E 
m c c 2 



= 2 - 5i °- ,s, jl) 2 Glv) w ' 



(8) 



where ctt is the Thompson scattering cross section, B the 
magnetic field strength and Ub its energy density. 

With the above equations, we can express the ratio 
of the diffuse synchrotron emission to the emission from 
SNRs as: 



-Psnr(^) 



N(E) dis (dE/dt) syn diS (dE/du) diS 
N(E) SNR (d£/di) synjSNR (d£/d^) SNR 
g^Wggfjff (d-E/d^) diff 

£-^TSNR£ 2 5 S 2 NR(d£/d^ SNR 



B, 



diff 



-BsNR 



Tdiff (^0) 
TSNR 



(9) 



Thus, the ratio of the diffuse synchrotron emission to the 
emission from SNRs is mainly determined by the ratio 
of the life-times of CREs in the two zones and by the 
ratio of the magnetic field strengths. The magnetic field 
strength in SNRs is expected to be higher than in the 
diffuse ISM because the gas in which the magnetic field 
is frozen-in gets compressed by the shock. In general, the 
magnetic field strength found in SNRs exceeds the value 
expected by this compression which would be a factor of 4 
in a strong shock, indicating that the magnetic field is not 
only compressed, but enhanced by the shock. Especially 
in young SNRs, the inferred values span a wide range from 



a few [iG (SN 1006; Tanimori et al. 1998) to much higher 



values of 70 fiG (Kepler's S NR; M atsui et al. |1984[ ), > 80 
jj,G (Cas A; Cowsik & Sarkar |l980| ) and even 300 fiG (Crab 
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Nebula; Kennel & Coroniti 1984). Taking these values as 
a guide, we assume for our present estimate Bsnr = 75/iG 
and Bdiff = 5/xG. Furthermore, we adopt for the injection 
spectral index 7 = 2.2 (Volk et al. 1988), for the life-time 
of a SNR its adiabatic phase tsnr ~ 3 10 4 yr and for r 
the energy loss time scale due to synchrotron and inverse 
Compton losses. The latter yields for a magnetic field of 
5 /xG and a radiation field of energy density 1 eV cm" 3 
Tdiff = 2.5 10 7 yr at v = 1.5 GHz. With these values we 
get -PsNR./-Pdiff = 0.09 at a frequency 1.5 GHz. 

A further contribution to the radio emission of a galaxy 



are radio supernovae (RSNe, e.g. Weiler et al. |1986[ ). Their 
contribution is difficult to estimate due to uncertainties in 



their luminosity and life-time. Perez-Olea & Colina (1995) 
have included RSNe in a model for the radio emission of 
starburst galaxies by generalizing and extrapolating the 
luminosities of 9 observed RSNe. Adopting their numbers, 
we calculate that for a galaxy in a steady state the total 
radio emission of RSNe is, for a life-time of a RSNe of 100 
yr, 1.5 x the total radio emission of SNRs. For a shorter 
life-time of 10 yr, their total radio emission would be a 
factor of 2 lower. Thus, in spite of the uncertainties, it 
seems plausible that RSNe contribute at a level of at least 
a few percent to the total radio emission of a galaxy. 

2.2. Comparison to observations 

In order to compare this theoretically derived ratio to ob- 
servations we calculate the expected average radio emis- 
sion of a SNR: 



= 1.0710" 



-25 / 9SN 



eV" 



-) WHztlO) 



MHz 



where we have inserted the relations from eqs. (|6|) and (g). 
For 7 = 2.2 this predicts a radio spectral index of SNRs of 
0.6, in reasonable agreement with the observations. The 
total number of CREs per energy interval produced by a 
SN, ggNj can be estimated from the total energy released 
by a SN which we take as 10 51 erg. Assuming that 10% of 
this energy goes into CR acceleration and of this 1% into 
the electron component, a total energy of Pg N = 10 48 erg 
is available per SN to produce the relativistic electrons 
following a power-law spectrum: 



^sn = <Zsn y E E {-) AE. 



(11) 



With 7 = 2.2, assuming for the lower energy limit of the 
produced spectrum E m i n = 1 GeV and adopting the scal- 
ing parameter Eq = 1 GeV this yields 



<?SN 



El 



En 



7-2 



(7 - 2) = 1.25 x 10 41 cV _1 . (12) 



In total, we obtain for the radio emission of a SNR, 
taking Psnr = 75/iG, 



Psnr = 2.1 10 1 



GHz 



WHz~ 



(13) 



This can be compared to the average observed emission 
of a SNR which can be derived from the relation between 
the surface brightness £ and the diamete r D o f SNRs. We 
adopt the relation found by Huang et al. ( 1994 ) who fitted 
SNRs from the Galaxy, M 82 and the Magellanic Clouds: 



-3.6 



I.v kui,. - 6.0 Mr 16 ( — J Wm^Hz^sr- 1 . (14) 



Following Condon & Yin (1990) one can connect the 
diameter and the age of the radio emitting SNR by: 



^U0.43 (1 

pcy V n J \y r 



(15) 



where P50 is the energy of a SN going into particle ac- 
celeration in units of 10 50 erg and n is the gas density in 
particles per cm 3 . In the following we assume P50 — 1 and 
n = 1. Using L = 47r 2 (P/2) 2 I] we then calculate the av- 
erage luminosity of a SNR by integrating the luminosity 

evolution from t = to t = Tsnr — 3 10 4 yr and adopting 

-0.6. 



the frequency dependence of eq. (rL3[) , Psnr oc v 



< Psnr >= 4tt 

c3 10'" 



/ 1 



V T SNR 
17 / v 



I s ^ = 3 - 31017 (gh^ 



WHz~ 



(16) 



which is very close to the above theoretical value (eq. 13). 

We can now estimate the total emission of SNRs ex- 
pected for our Gal axy: A dopting a SN rate of 1 SN per 
30 yr (Berkhuijsen 1984 ), the total number of the SNRs 
present in our Galaxy is, with tsnr = 3 10 4 yr, A^snr = 
^SN^sn = 1000. This yields, together with eq. ([IJ]), a 
total radio emission from SNRs of PsNR,tot(l-5 GHz) = 
P SN r(1.5GHz)A snr = 2.6 lO^WHz" 1 . This can be com- 
pared to the total radio emission of our Galaxy which is 
at 408 MHz 8.7 lO^WHz" 1 (Beuermann et al. |1985[ ). Ex- 
trapolating this value, with a radio spectral index of 0.8, 
to 1.5 GHz we get 3 10 21 W Hz -1 . From this calculation 
we derive a contribution of about 9% from SNRs to the 
total radio emission, in satisfactory agreement with the 
value derived in the previous section in spite of the crude- 
ness of both estimates. Therefore, we conclude that ~ 10% 
is a reasonable estimate for the contribution of SNRs to 
the total radio emission. This value is also in agreement 
with previous estimates by other authors, who compared 
the expected total synchrotron emission from SNRs with 
the observed radio emission of galaxies, and have found 
that SNRs alone can only expla in about 10% o f the total 
radio emission (Biermann 1976 ; Ulvestad 1982 ). 
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Fig. 1. The nonthermal spectral index, a n thj as a function 
of the diffuse radio spectral index, adiff, between 1.5 and 
5 GHz, for different contributions of the radio emission 
from SNRs to the radio emission: /snr =5% (full line), 
10% (dashed line), 20% (dotted line). 



2.3. Influence on the spatially integrated radio spectral 
index 

In Fig. [| we show the total (i.e. spatially integrated) non- 
thermal spectral index, a n th, as a function of the spectral 
index of the diffuse radio emission, adiff> f° r different con- 
tributions of SNRs to the radio emission, expressed by 
/snr = -PsNR/fnth- We assume, as above, an (energy) in- 
jection spectral index for CREs of 7 = 2.2, yielding a (fre- 
quency) spectral index of SNRs of cusnr = (7 — 1)/2 = 0.6 
(eq. (|10|)). A lower source spectral index of 7 = 2.0 to 2.1 
i s ind icated by recent calculations by Berezhko and Volk 
( |l997| ). This would give asNR = 0.5 to 0.55 and would 
somewhat increase the difference between a n th and adig 
for a given /snr- 

The radio emission of SNRs has a noticeable effect on 
the nonthermal radio spectral index. Due to the contri- 
bution of SNRs it gets natter, especially for large values 
of adiff where the spectral index is lowered by ~ 0.1 for 
/snr =0.1. 



3. CR propagation and energy losses 

The types of propagation of CREs in a galaxy are dif- 
fusion, i.e. the CREs are scattered randomly by magnetic 
field irregularities, and convection, a systematic movement 
of the scattering centres outwardly. While propagating 
through a galaxy, CREs suffer energy losses. Eventually, 
CREs might escape from the galaxy. 



In order to discuss the radio spectral index in the var- 
ious cases, we calculate two simple models, one that con- 
siders diffusion and one in which we study convection. In 
both models we take into account energy losses by syn- 
chrotron and inverse Compton losses. Here, we ignore adi- 
abatic and bremsstrahlung losses which might also be im- 
portant, especially in starburst galaxies, but since they do 
not change the spectral index, they are in the framework 
of the present study equivalent to no energy losses at all. 

3.1. Diffusion model 

The diffusion model can be described by the following 
equation: 



D{E) 



d 2 N(E, z) d 



+ — (b(E)N(E,z)) = -Q(E,z), (17) 



dz 2 ' dE 

where z is the height above the galactic midplane. We as- 
sume that the sources of CRs are situated in the midplane 
at z — and that their energy dependence is given by a 
power-law according to eq. (|2|). The synchrotron and in- 
verse Compton losses, b(E), are assumed to be spatially 
homogeneous and are given by: 



b(E) 



dE 



-arc 



syn+iC 



E 



m c c 



rad 



(18) 



t/rad is the energy density of the radiation field below the 



Klein- Nishina limit (e.g. Longair 1992). We adopt a free 
escape boundary at \z\ — Zh, the outer boundary of the 
halo, so that N(E, \zh\) = 0. An approximate solution is 
given by (for further details see Lisenfeld et al. 1996a): 



N(E, z) = v m 



E 



m e c 2 J \lGeV/ 



iF 1 (p+ | i | 1 Sh 
1-F1O, 5,s h ) 



- iFi(p+^,s)-\z\ 



with 



P 



2x + /1 — 3 



(19) 



(20) 



2(1 -m) J ' 

iFi(a,6, x) denoting the confluent hypergeometric func 
tion (or Kummer function), and 

Cl(Efrad + Efc)(l-/*) 2 ( E \ ( E \ 



4Dn 



m c c' 



2 ) VlGeVy 



(21) 



and 



4.Dn 



m e c 2 / \ lGeV 



where 
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4 cttc 



3 m e c 2 



(23) 
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3.2. Convection model 

The convection model is described by the following equa- 
tion: 



1.1 



Q(E,z) = 



§- z (VN(E,z)) 



_d_ 

dE 



18V 
3~dz~ 



E + b(E) \ N(E,z) 



(24) 



where V is the convection speed, assumed here to be con- 
stant. The solution of this equation can be obtained with 
the method of characteristics: 



N(E,z) = -1-E-i 



1 



bEz 



(25) 



3.3. Asymptotic spectral indices 



The diffuse synchrotron emission is calculated from 
N(E, z) in the same way as for SNRs (eqs. (0) and (§)). 
The total diffuse synchrotron emission is obtained by inte- 
grating Pdiff over the extent of the halo, Zhaio, which in the 
case of the diffusion model coincides with the free-escape 
boundary, z^. 

We define the escape probability, P eS c(v)i as the ratio 
of the measured total nonthermal radio emission (includ- 
ing the diffuse synchrotron emission and the radio emis- 
sion from SNRs), P n th(v), to the maximum possible non- 
thermal radio emission, ■P n th,max( l/ )) that would be mea- 
sured if the halo were infinitely large. 

The expected diffuse and nonthermal spectral indices 
for different asymptotic cases are: 

Case (1): Synchrotron and inverse Compton losses are 
very strong so that no escape is possible. Then, indepen- 
dent of whether diffusion or convection is the principle 
mode of propagation, we have adiff = 7/2 = 1.1. Includ- 
ing the contribution of SNRs, the spectral index becomes 
(see Fig. @) a nt h — 1.0. For a lower injection spectral in- 
dex, 7 = 2.0, we would get a n th — 0.9. 

Case (2): CREs can escape freely from the galaxy by 
diffusion without suffering synchrotron and inverse Comp- 
ton losses. The diffusion coefficient can be energy depen- 
dent, D(E) oc E^. For our Galaxy, /j, has been determined 
from the energy-dependence of the fraction of secondary 
to primary CR nuclei. For CRs with energies above a few 
GeV, the energy range relevant for radio observations in 
the GHz range, [i = 0.4 — 0.7 has been found (Garcfa- 
Muhoz et al. 1987). The expected diffuse spectral index 
is then, adopting [i = 0.5, adiff = (7 + M — l)/2 = 0.85, 
and including the contribution of SNRs it becomes a n th — 
0.82, with corresponding reductions for 7 = 2.0. 

Case (3): CREs propagate by convection and do not 
suffer considerable synchrotron or inverse Compton losses. 
In this case we expect adiff = ctnth = (7 — l)/2 = 0.6 
for 7 = 2.2; the spectral index remains unchanged when 
including the emission from SNRs. 
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Fig. 2. The nonthermal spectral index between 1.5 and 
5 GHz, assuming 7 = 2.2 and /snr = 0.1, is shown as 
a function of the escape probability P sc- The full line is 
calculated for a convection model, the dashed line for a 
diffusion model with /1 = 0.2 and the dotted line for a 
diffusion model with /1 = 0.5. The full and dashed-dotted 
horizontal lines give the observational value for a n th and 
its dispersion, calculated as the average of the results by 
Niklas et al. (gggg) and Klein (fljggl). For 7 = 2.0 all model 



curves would be lower by about 0.1. 



3.4- The spatially integrated spectral index 

The results of the models are shown in Fig. where the 
nonthermal spectral index (assuming a contribution of 
SNRs of 10% at 1.5 GHz) is plotted as a function of the 
escape probability for the case of diffusion and of convec- 
tion. The values for a n th predicted by the two models are 
compared to the observed mean nonthermal spectral in- 
dex which is about 0.85 ± 0.04. This observational value is 
obtain ed as the average of the results fr om re sults by Klein 
( |1988D (0.88 ± 0.06) and Niklas et al. ([L997D (0.83 ± 0.02) 
who determined the average a n th of samples of galaxies by 
subtracting the thermal radio emission from the observed 
total radio emission. 

Within the pure diffusion model with /1 = 0.5 the 
observed a n th is explained by a high escape probability 
(P csc > 50%). For a weaker energy dependence of the dif- 
fusion coefficient, /1 = 0.2, however, the predicted escape 
probability is lower, P csc = 30 — 60%. The conclusion from 
the convection model is very different: Here the observed 
a n th is explained by P csc ~ 25%. The different conclusions 
with respect to the escape probability given by the diffu- 
sion and the convection model are due to the fact that 
(i) convection is an energy independent process and (ii) 
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Fig. 3. The diffuse radio spectral index as a function of the 
height over the galactic plane, z, in the case of important 
energy losses. The full line is calculated for a convection 
model, the dashed line for a diffusion model with /i = 0.5 
and the dotted line for a diffusion model with n = 0.2. We 
have assumed a magnetic field B = 5 /iG, a radiation field 
of energy density 1 eV cm~ 3 , convection velocity V = 200 
km s _1 and diffusion coefficient Dq = 10 29 cm 2 s _1 . 

in the case of convection the relation between the spectral 
index and the escape probability is not linear, in the sense 
that low escape probabilities are predicted for almost the 
whole range of spectral indices. 

We conclude that from the observed average radio 
spectral index of galaxies it can neither be decided to 
which degree CREs can escape from galaxies nor whether 
diffusion or convection is the dominant mode of propa- 
gation. Only for asymptotic cases partial conclusions can 
be drawn: If the spectral index of a galaxy is very steep 

0.9) energy losses are important and the escape rate is 
low. Yet, it cannot be determined whether diffusion or con- 
vection is the dominant type of propagation. A very flat 
overall spectral index (< 0.7), on the other hand, indicates 
strong convection or almost energy independent diffusion. 
Since measurements in our Galaxy indicate that diffusion 
does depend significantly on energy (fi = 0.4 — 0.7, Garcia- 
Muiioz et al. 1987), convection is the more likely expla- 
nation. In this case, however, we cannot draw any conclu- 
sions with respect to the confinement of CREs. 



3.5. Spatially resolved spectral index 

The question of how important energy losses are for the 
propagation of CREs in a galaxy can be decided if the 
spectral index distribution in the halo is known. If the 
spectrum steepens considerably away from the disk, it is 
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Fig. 4. The diffuse radio spectral index as a function of the 
height over the galactic plane, z, in the case of almost free 
escape. The line coding and parameters are as in Fig. ||. 



a clear sign that CREs suffer considerable synchrotron and 
inverse Compton losses as they propagate outwardly. This 
is illustrated in Fig. || and |J, where we show the diffuse ra- 
dio spectral index, adiff, z , calculated from the convection 
and diffusion models of Section 3.1, as a function of the 
distance from the galactic plane, z. In Fig. || we show the 
results for a large halo (z^ = 15 kpc) in which the energy 
losses affect the CREs considerably before they manage 
to escape. The escape probabilities are here 9% for con- 
vection and 6% (fi — 0.2), respectively 10% = 0.5) for 
diffusion. In Fig. |]thc opposite case is shown: A small halo 
(z/j = 3 kpc) where the CREs can escape without suffering 
significant energy losses, so that the escape probabilities 
lie between 78% for convection and 85% (fi — 0.2), respec- 
tively 87% (/i = 0.5) for diffusion. We find that important 
energy losses are reflected clearly in a steepening of the 
spectral index with increasing distance from the plane. 
On the contrary, if the CREs can escape almost freely, or 
if only adiabatic or bremsstrahlung losses occur, no such 
steepening is seen. 

3.6. Comparison to observations: Are energy losses 
important? 

The distribution of the spectral index in the halo is a 
definite diagnostic of whether energy losses are impor- 
tant. The observational study of radio halos in external 
galaxies is difficult, because the halos are extended and 



intrinsically faint (see Dahlem 1997). Therefore, not many 
galaxies with multi-frequency radio data for the halo are 
known. 
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Fig. 5. A comparison of our model results to the data of 



NGC 891 (Hummel et al. \L99l\ ). The triangle refer to the 
spectral index in the western part of the halo, and the stars 
to the eastern side of the halo. The full line corresponds to 
the convection model and the other lines to the diffusion 
model with /i = 0.5 and different values for the diffusion 
coefficient and thus escape probability (see text). 



One of the best studied radio halos is that of the spiral 
galaxy NGC 891 (Hummel et al. |l99l|) . Its spatially intc- 
gr ated n onthermal radio spectral index is 0.78 (Niklas et 
al. 



1997 ), only slightly below the average spectral index of 
galaxies (see Fig ||). In Fig. || we show its radio spectral 
index distribution in the halo. The data are at a resolution 
of 40" , corresponding to 1.8 kpc at the distance of 9.5 kpc 
of NGC 891. The spectrum steepens very quickly within 
the first 2-3 kpc. It continues to steepen more slowly out- 
side this range, especially in the eastern side of the halo. 
Here, the steepening becomes more pronounced at z Z 6 
kpc. 

In Fig. H we also show the results of our model in 
which we have included, in order to correspond to the 
data, apart from the diffuse radio emission, the contribu- 
tion from SNRs and the thermal radio emission. For the 
latter two components we have assumed, /snr = 0.1, as 
before, and a contribution of the thermal radio emission of 
10 % at 1.5 GHz, a typical value for normal galaxies (Con- 
don 1992). We furthermore assume that the thermal radio 
emission and the radio emission from SNRs is restricted 
within 100 pc above the disk (the exact value of this spa- 
tial extent is not important because of the low resolution 
of the data). Finally, we convolve the model results with 
a Gaussian beam of HPBW 1.8 kpc, corresponding to the 
resolution of the data. 



We have adjusted the diffusion coefficient and convec- 
tion velocity such that we achieved the best possible fit 
to the data. Because of the simplicity of the models de- 
scribed here, we cannot expect to describe the data in all 
detail, we are aiming rather at a qualitative comparison 
with the main features. For a quantitative comparison a 
more complex model would be necessary, combining dif- 
fusion and convection, taking into account spatial varia- 
tions of the parameters (B,D , etc.) and describing more 
realistically the distribution of SNRs and the contribution 
of the thermal radio emission. Nevertheless, the diffusion 
model is able to describe a large range of the data reason- 
ably well. The fast steepening of the spectral index within 
the first 2-3 kpc above the disk can to a large extent be 
attributed to the contribution of SNRs and the thermal ra- 
dio emission close to the disk. The steepening outside this 
range requires the presence of dominant energy losses: The 
diffusion models that describe the data well have escape 
probabilities of 11 % (dotted line) and % (dashed line). 
For comparison a model with a high escape rate is shown 
(78 % - dashed-dotted line) which is not able to describe 
the data. In order to reproduce the fast steepening of the 
spectral index at z > 6 kpc in the western part of the halo, 
we would possibly have combine the diffusion model with 
the convection model. 

Radio halos have been observed in several other galax- 
ies, although mostly not in such detail. In all of these 
galaxies a steepening of the spectral index at least up 
to ctnth ~ 1.0 is found with increasing distance from the 
disk. This is the case for both starburst galaxies, as M 82 



(Seaquist & O dega rd |199J Reuter et al. |1992[ ); N GC 25 3 
(Carilli et al. |l992|) ; NGC 2146 (Lisenfeld et al. |l996b|) ; 
NGC 4666 (Sukumar et al. |l988| ; Dahlem et al. |l997| ) and 
non-starburst spiral galaxies like NGC 4631 (Hummel & 
Dettmar |1990|) ; NGC 5775 (Duric et al. |l998| ); NGC 5055 
and NGC 7331 (Hummel & Bosma |l982[ ). We conclude 
therefore that for all galaxies for which a multi-frequency 
radio halo has been observed, energy losses are important 
and it seems likely that this is generally the case in galactic 
halos. 



Niklas & Beck (1997) have found a trend that actively 
star-forming galaxies tend to have lower spectral indices 
and more quiescent galaxies steeper ones. They interpret 
this as an indication that in galaxies with a high star- 
formation efficiency (SFE) CREs escape more easily due 
to a galactic wind. On the basis of the above considera- 
tion, however, it seems more likely that in these galaxies 
convection causes the flat spatially integrated spectrum. 
Energy losses play an important role also in galaxies with a 
high SFE which is shown by the steepening of the spectral 
index in the halo. A good example is the starburst galaxy 
M 82. For this galaxy the observed overall synchrotron 
spectral index a n th = 0.66 (Niklas et al. 1997) indicates 
strong convection (see Fig. 2) whereas at the same time 
the steepening of the spectral index in the halo shows the 
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presence of important energy losses (Seaquist & Odegard 



1991; Reuter et al. 1992) 



3.7. Consequences for the FIR/radio correlation 

The question whether CREs can escape freely from galax- 
ies or whether synchrotron and inverse Compton losses 
determine their spectra has been asked in the framework 
of the FIR/radio correlation. Volk ( 1989) ) has proposed a 
calorimeter model in which the tightness of the correla- 
tion can be explained by the fact that CREs loose their 
energy by synchrotron and inverse Compton losses below 
the energy level corresponding to the observing frequency 
before they a re able to escape from a galaxy. Lisenfeld 
et al. ( 1996a ) have generalized this model, allowing for 
moderate escape and a spatially inhomogcncous magnetic 
field. Other authors have claimed the opposite, i.e. that 
CREs can esc ape m ore or less freely from the galaxy (Chi 
fc Wolfendale |l990t Helou & Bicay |l993|; Niklas & Beck 



19971 ) 



It has been argued (Niklas & Beck 1997 ) that the ob- 
served spectral index of galaxies, a n th — 0.85, is in contra- 
diction with the calorimeter model which predicts in its 
asymptotic version a n th — 1.0—1.1, corresponding to case 
(1) of Section 3.3. The difference between the prediction 
of the calorimeter model and observations of the radio 
spectral index can be decreased by allowing for moder- 
ate escape and a spatially inhomogeneous magnetic field 



(Lisenfeld et al. 1996a) and even more by including the 
contribution of SNRs (this work), so that such a modified 
calorimeter model would predict a nt h — 0.85 — 0.95, cor- 
responding to the observed average value. However, even 
in this case this model is not able to account for the flat- 
ter spectra that some of the galaxies (e.g. M 82) show 
because it is based on diffusion. In this paper, we have 
shown, that for galaxies with such a low spectral index, 
convection is likely to be important. With respect to the 
question of whether confinement or escape is taking place 
the overall spectral index is not a good diagnostic. Spa- 
tially resolved observations are necessary and these indi- 
cate that energy losses indeed do play a dominant role. 
Thus, even for galaxies with a low overall spectral index, 
like M 82, escape seems to be negligible as predicted by 
the calorimeter model. 

4. Conclusions and summary 

We have estimated and discussed different processes and 
components that shape the radio spectral index of a 
galaxy. In particular, we have considered the contribution 
of SNRs, the influence of different types of CR propaga- 
tion (diffusion and convection) and of energy losses (in- 
verse Compton and synchrotron) . Our conclusions can be 
summarized as follows: 

(i) We have estimated that the radio emission of SNRs 
represents about 10% of the nonthermal radio emission 



of a galaxy. This moderate contribution has a noticeable 
effect on the nonthermal radio spectral index, lowering it 
by 0.1 for steep spectra. 

(ii) The spatially integrated radio spectral index of a 
galaxy is difficult to interpret, as it is influenced by various 
processes. Only for asymptotic cases conclusions can be 
drawn: In galaxies with very steep spectral indices (^ 0.9) 
energy losses are high and the escape rate of CRE is low. 
A very flat spectral index (< 0.7), on the other hand, is a 
sign for strong convection. 

(iii) In order to draw more general conclusions about 
the CR propagation and the importance of energy losses, 
spatially resolved radio observations of galactic halos are 
necessary. A steepening of the spectral index away from 
the galactic disk is a clear indication that inverse Compton 
and synchrotron losses are important. 

(iv) In all galaxies with multi-frequency radio data for 
the halo such a steepening has been found and therefore 
in these galaxies synchrotron and inverse Compton losses 
take place. Thus, it seems likely that these energy losses 
are generally important in galactic halos and that escape 
rates are low. The low overall spectral index that is fre- 
quently found in star burst galaxies, at the same time as 
the steepening of the spectrum away from the disk, can 
be interpreted as convection. 
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